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Abstract 

We present the complete NLO electroweak contribution of 0{c? s oi) to the production of diagonal 
squark-anti-squark pairs in proton-proton collisions. Compared to the lowest-order 0(a s a + a 2 ) 
electroweak terms, the NLO contributions are also significant. We discuss the LO and NLO elec- 
troweak effects in cross sections and distributions at the LHC for the production of squarks different 
from top squarks, in various supersymmetric benchmark scenarios. LO and NLO can add up to 10% 
in cross sections and 25% in distributions. 
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1 Introduction 

The exploration of electroweak symmetry breaking will be one of the main tasks at the Large Hadron Col- 
, lidcr (LHC). Experiments are expected to either verify the Higgs mechanism of the Standard Model (SM), 

' or to detect signals of physics beyond the SM. The concept of supersymmetry [1,2] provides a promising 

alternative version of the Higgs mechanism where symmetry breaking occurs without introducing new 
scalar couplings that potentially can become strong, thus stabilizing the electroweak scale. The real- 
ization of supersymmetry in terms of the Minimal Supersymmetric Standard Model (MSSM) [3-5], has 
come up as the most promising extension of the SM, a predictive framework that allows to make precise 
predictions to be investigated by indirect and direct experimental studies. The indirect access through 
virtual effects in electroweak precision data [6] provides an overall fit [7,8] with a quality at least as good 
as in the SM, in specific observables like g — 2 of the muon [9, 10] even better, and yields bounds on the 
light Higgs boson mass with less tension than in the SM [11-13]. 
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If supersymmetry (SUSY) is realized at the TeV scale or below, it will be accessible to direct experi- 
mental studies at the LHC through the production of SUSY particles. In particular, colored particles like 
squarks and gluinos will be copiously produced, and the hadronic production of squark- anti-squark pairs 
is expected to play an important role for SUSY hunting. The cross section is in the range from 0.5 to 10 
pb for masses of squarks and gluinos below 1 TeV and can be measured with a statistical uncertainty of 
the order of a few percent even in the low luminosity regime. Moreover, squark cascade decays into q\\ 
lead to a signature with missing Et plus jets and possibly leptons that is well suited to detect MSSM 
signals [14, 15]. The number of hard jets allows the distinction between gluino and squark decays. Finally 
with the help of decay chains one can reconstruct the mass of the squarks up to 2 TeV with a resolution 
better than 10 % [16-18]. 

The first prediction of the cross section for hadronic production of squark pairs in the early 1980's 
was done at lowest order 0(a 2 s ) in supersymmetric QCD [19-23]. QCD contributions at NLO, 0(a^), for 
the processes PP — > Q a Q b *X, PP — ► Q a Q b X (Q ^ t) were calculated more than ten years later [24,25]. 
They increase the cross section by typically 20 to 30 %, and they substantially reduce the dependence 
on the factorization and renormalization scale. NLO QCD corrections to the production of top-squark 
pairs, performed in [26], are also positive and can increase the cross section by 40-50%. 
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Besides the QCD-based production mechanisms, there are also partonic processes of electroweak ori- 
gin, like diagonal and non-diagonal squark pair production from qq annihilation [27,28]. They proceed 
through s-channel photon and Z exchange, and also through ncutralino/chargino exchange in the t- 
channel (if Q is different from t), yielding terms of 0(a 2 ) and O(a s o). Due to interference between the 
tree-level QCD and electroweak amplitudes for Q ^ t, the electroweak contributions can also become 
sizable, reaching values up to 20% [27]. 

For reliable predictions, electroweak contributions at NLO have to be taken into account as well. 
In the case of top-squark pair production [29-31], they were found to be significant, with effects up to 
20%. In general, NLO electroweak (EW) contributions consist of loop contributions to the tree-level 
amplitudes for qq annihilation and gluon fusion, together with real photon and gluon bremsstrahlung 
processes, yielding an involved struture of interference terms in qq annihilation. Moreover, photon-gluon 
induced parton processes also contribute owing to the non-zero photon distribution in the proton. In 
this paper we present the NLO electroweak contributions, of O(a 2 s o) 1 to the production of diagonal 
squark-anti-squark pairs different from top- and bottom-squarks, 

P P -> Q a Q a * X (Q ^ i, b) . (1) 

They show significant differences to top-squark production, based on the following pecularities. 

• In leading order 0(a 2 ), the squark pair Q a Q a * can be produced via annhilation of a Q Q pair 
through amplitudes that involve also the exchange of a gluino in the i-channel, thus enhancing the 
relative weight of the annihilation channel in (1). 

• Electroweak tree diagrams with t-channel neutralino and chargino exchange are part of the ampli- 
tudes for Q Q — > Q a Q a * and Q' Q' — > Q a Q a *, where Q' is the isospin partner of Q in a quark 
doublet, yielding EW-QCD interference already at the tree-level. 

• At O{a 2 s o) many types of interferences occur between amplitudes of O(a s o) and 0(a s ) as well as 
between 0(a 2 ) and 0(a) amplitudes. 

These features make the calculation of the EW contributions of 0(a 2 a) to the processes (1) more involved 
than in the case of it* production where no i-channel diagrams occur at lowest order. Our analysis shows 
that the EW effects of NLO can reach the same size as the tree- level EW contributions of O(a s o) and 
0(a 2 ), which we will include in our discussion as well. 

The case of bb* production will not be treated here. Owing to 6-tagging, bottom-squarks can be ex- 
perimentally distinguished from the squarks of the first two generations [18,32,33]. Moreover, in the case 
of bb production the partonic process bb — > bb exhibits specific features, like mixing between left- and 
right-handed 6-squarks, mixing angle rcnormalization [34] , non- negligible Higgs-boson contributions and 
enhanced Yukawa couplings for large vaues of tan/3 with the need of resummation [35]; other peculiarities 
for massive initial-state partons are the proper counting of the orders of the perturbative expansion [36,37] 
and the appropriate choice of the factorization scale [38] . A dedicated extra analysis for 6-squark final 
states thus seems appropriate. 

The outline of the paper is as follows. In Section 2 we briefly summarize the various tree- level 
contributions to the processes (1). Section 3 describes the structure of the NLO terms of EW origin that 
contribute at 0{a 2 a) and the strategy of the calculation. Evaluation of the EW effects and their analysis 
for the LHC are presented in Section 4 and summarized in Section 5. A list of Feynman diagrams and 
counter terms with specification of rcnormalization, and technical details for the calculation of singular 
integrals are collected in the Appendix. 
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2 Tree-level contributions to squark pair production 



In this section we list the lowest-order cross sections for the process (1) arising from tree-level amplitudes 
at order 0(a 2 ), 0(a s a) and 0(a 2 ). We will use the convention da^ to denote the cross section for 
a partonic process X at a given order 0(a"a b ) in the strong and electroweak coupling constants. The 
parton luminosities for getting to the hadronic cross section are given by the convolution 

dLij , „ 1 f 1 dx 



dr 1 + 5y 

where fi(x) is the momentum distribution of the parton i in the proton (PDF). 
2.1 Squark pair production at leading order 

The leading-order contribution to the process (1) is QCD based, of 0(a 2 ). In the notation mentioned 
above, the differential cross section reads as follows, 



(3) 



The sum runs over the quarks q = u, d, c, s. S and s = tS are the squared CM energies of the hadronic 
process (1) and of the partonic subprocess, respectively. Moreover, with the squark mass ttIq q , the 
threshold value t is determined by r = 4mj~ JS. 

da 2 ! ~ - and da 2,0 ~ - denote the 0(a 2 ) differential cross sections for the partonic processes 

qq^Q a Q a * gg^Q a Q a * y s/ 

q(pi)q(P2) - Q a (h) Q Q *(fc 2 ), (4) 
g(pi)g(P2) - Q a (ki) Q a *(k 2 ), (5) 

respectively, which are obtained from the Feynman diagrams in Fig. 1 of Appendix A. Explicit expressions 
for these leading-order cross sections can be found in Refs. [19,22,23]. Owing to flavour conservation in 
SUSY QCD, the diagram with the exchange of a gluino in the t channel contributes only if q = Q. 

2.2 Tree-level electroweak contributions of 0{ol s ol) and 0(ct 2 ) 

The 0(a s a) and 0(a 2 ) contributions to the process (1), involving electroweak terms, can be written as 
follows, 

q t o ^ j 
2 

The parton cross section da q L^Q a Q a , is obtained squaring the tree-level electroweak diagrams depicted in 
Fig. 2 of Appendix A. The diagram with i-channel neutralino exchange contributes only if q = Q, and the 
diagram with chargino exchange appears only if q' = Q, q' being the SU(2) partner of the quark q, since 
we treat the CKM matrix as unity, da 1 '} - ~ originates from interference between the aforementioned 
tree-level electroweak diagrams and the tree-level QCD graphs of Fig. 1. Analytical expressions for these 
cross sections can be found in Ref. [27]. 

As a new element at 0(a s a), photon-gluon fusion occurs as a further partonic process, 

7(Pi) fl(P2) -> Q a (ki) Q a *(k 2 ). (7) 
The corresponding cross section, with t = (pi — fci) 2 , 
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contains the spin- and color-averaged squared tree-amplitudes from the diagrams in Fig. A, 



— — m% (m% + s 2 ) + tu(tu — 2m% ) 

i i >0 |2 -,c 2 Q,a y Q.a ' v Q,a' 

1 7S^Q a Q" 1 S Q {t-m% ) 2 (u-m 2 ~ ) 2 

with u = (pi — fc 2 ) 2 , and the electric charge of the squark Q. The presence of photons in the proton 
follows from including NLO QED effects in the evolution equations for the PDFs. The photon PDF is 
part of the publicly available PDF set of [39]; together with the gluon PDF, the 717 luminosity entering (6) 
is built according to Eq. (2). 

3 Virtual and real 0(ot 2 s a.) corrections 

In this section we describe the computation of the 0{a 2 a) corrections to the process (1) arising from 
loops and from photon/gluon bresmsstrahlung. The corresponding contributions to the hadronic cross 
section are expressed in obvious notation, 

, cw,NLO f 1 1 dLgg / 2,1 , J 2.1 \ 




2,1 



dT 

T 



dL„„ ,.,51 dLn, 



qg ( \ , 2,1 . aJ ^qg 1 \ 1 2,1 



(9) 



dr qg^Q a Q a 'q dr qg^Q a Q a 'q 

Other bremsstrahlung contributions to the hadronic cross section are of the type 

7(Pi) 9(P2) - Q a (ki) Q a *(k 2 ) q(k 3 ), 7 (pi) KV2) - Q a (h) Q a *(k 2 ) q(k 3 ). (10) 

We will not consider this class of processes here; they are further suppressed by an additional factor a s 
with respect to process (7), and thus negligible. 

Diagrams and corresponding amplitudes are generated using FeynArts [40,41]. The algebraic treat- 
ment and numerical evaluation of loop integrals is performed with support of FormCalc and LoopTools [42, 
43]. IR singularities are regularized by a small photon mass, while quark masses are kept as regulators 
for the collinear singularities. 



3.1 Gluon fusion with electroweak loops 

The first class of corrections entering Eq. (9) are the O(a 2 o) electroweak virtual contributions to gg 
fusion (5), given by the partonic cross section 

A4° is the tree level gg amplitude (Fig. 1), and .M 1 '*™ is the one-loop amplitude with EW insertions in 
the QCD-based gg tree diagrams. These loop diagrams do not depend on the flavour of the final squark 
and thus they are identical to those listed in [30] for the particular case of Q — i. We therefore do not 
repeat them here. 

In order to get rid of the UV divergences we have to include the proper countcrterms for one-loop 
renormalization. Their explicit expressions in terms of the renormalization constants can be found in 
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Appendix B. In the case of the gluon- fusion subprocess we have to renormalize the squark sector only. 
We use the on-shell scheme [30,34,44], where the independent parameters for a squark isospin doublet 1 
are chosen to be the masses of the two up-squarks, the mass of one of the two down-squarks, and the two 
mixing angles (which, however are irrelevant for the light-quark squarks where mixing can be neglected) . 
The actual expressions for the renormalization constants are also given in Appendix B. 

Notice that part of the virtual corrections to squark pair production are loop diagrams for the gluon- 
gluon-7J° vertex, with the heavy neutral MSSM Higgs boson H°. These terms become resonant when 
m H a > 2rriQ a and have to be considered a contribution to the process of H° production via gluon fusion 

with the subsequent decay H° — > Q a Q a *, rather than an electroweak loop correction. We will not consider 
scenarios in which such resonances occur. 

3.2 Gluon fusion with real photon emission 

The IR singularities arising from virtual photons in (11) are cancelled by including bremsstrahlung of 
real photons at 0(ala), 

g(pi) g(P2) - Q a (h) Q a *{k 2 ) 7 (fc 3 ) , (12) 

according to the diagrams depicted in Fig. 4. The integral over the photon phase space is IR divergent 
in the soft-photon region, i.e. for fc§ — > 0, and cancels the corresponding virtual singularities when added 
to the virtual contributions according to Eq. (9). 

For the technical treatment of photon-momentum integration and isolation of divergences we apply 
two different procedures: the methods of dipole subtraction and of phase space slicing. In the dipolc 
subtraction approach, one has to add and subtract an auxiliary function to the differential cross section 
that matches pointwise the singuilarities and is easy enough to be integrated analytically; the integral 
over the subtracted cross section is convergent and can be done numerically. Due to the universality of the 
soft singularities general expression for these functions are available. In particular we use the expressions 
given in Ref. [45]. Although the formulae quoted in this reference apply to processes involving fermions 
only, they can be generalized to processes with charged bosons owing to the universal structure of the IR 
singularities. 

The phase space slicing technique restricts the phase space integration to the region with a minimum 
photon energy AE = S s ^/s/2. The integration over this region is thus convergent and can be performed 
numerically. The complementary integral over the singular region with k% < AE can be done analytically 
in the eikonal approximation [46] , which is a good approximation if the cut S s is sufficiently small. More 
details are given in Appendix C. Comparison between the two methods provides a non trivial check of 
the computation. As illustrated in Fig. 12, the two methods yield results which are in good numerical 
agreement. 

3.3 qq annihilation with electroweak and QCD loops 

The structure of the parton processes of qq annihilation at higher order is more involved and requires a 
simultaneous treatment of electroweak and QCD loops. The virtual contributions of one-loop order to 
the partonic cross section is given by the interference of tree-level and loop amplitudes, 

where A / f°' qcd (M a '°™) is the amplitude related to the tree-level QCD (EW) diagrams depicted in Fig. 1 (2). 
A^ 1 '™ is the one- loop amplitude arising from the EW corrections to the QCD tree- level diagrams and the 

x Due to SU(2) invariance, mass renormalization of the different squarks from the same SU(2) doublet is correlated and 
has to be performed simultaneously. 
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QCD corrections to the EW tree-level diagrams. Finally, .A4 1 ' qcd is the one-loop amplitude corresponding 
to the QCD corrections to the QCD tree-level diagrams. 

The diagrams entering M 1 '"™ are displayed in Figs. 5-7 of Appendix A. They also contain the dia- 
grams with counterterm insertions required for renormalization and cancellation of UV divergences. The 
counterterms and the necessary renormalization constants can be found in Appendix B. Besides squark 
renormalization, also quark renormalization is needed. 

A4 1,qcd can kg obtained from the Feynman diagrams in Fig. 8 of Appendix A, including the proper 
counterterms. Besides renormalization of the squark sector, we have to renormalize also the gluino mass, 
the strong coupling g Sl and the quark-squark-gluino coupling g Sl which is related to g s via supersymmctry. 
The strong coupling constant is renormalized in the MS scheme, modified according to Ref. [24] in order 
to decouple heavy particles (top, gluino, squarks) from the running of a s . For the non-standard loop 
contributions, this procedure is equivalent to the "zero- momentum subtraction" used in Ref. [47]. Since 
dimensional rcgularization violates supersymmctry in higher orders, a finite difference between g s and g s 
is encountered at one-loop order. Supersymmetry is restored by shifting the renormalization constant for 
g s by the corresponding finite amount, which means an unsymmetric renormalization of g s and g s . More 
details and the specification of the counterterms can be found in Appendix B. 

3.4 qq annihilation with real photon emission 

The diagrams in Fig. 9 of Appendix A constitute the generic amplitude for photon bremsstrahlung at 
0(a 2 a) in the qq annihilation channel, 

q(Pi) q(P2) - Q a (fci) Q a *(k 2 ) 7 (fc 3 ) • (14) 

The corresponding cross section is singular both in the IR soft-photon region and in the collinear region 
(e.g. whenever k 3 pi — ► 0). Although IR singularities cancel in sufficiently inclusive observables, collinear 
singularities from initial-state radiation remain and have to be absorbed via factorization in the PDFs. 

The extraction of the singularities has been performed using the two different methods described in 
section 3.2. In phase space slicing, in this case, we have to introduce a further collinear cutoff 5 C on the 
angle between the photon and the radiating quark/antiquark. For sufficiently small S c , the integral over 
the singular region can be performed analytically. Explicit expressions can be found in Appendix C. In 
Fig. 12 we visualize the comparison between the two methods in the specific case of the partonic process 
uu — > u L u L *^/ as an example. 

3.5 qq annihilation with real gluon emission 

Finally, we have to take into account the class of qq annihilation processes with real gluon bremsstrahlung, 

q(pi) 5(P2) - Q a (h) Q a *(k 2 ) g(h) , (15) 

from either EW-based (Fig. 10 a) or QCD-based Born diagrams (Fig. 10 b). This class contributes to the 
cross section at 0(a 2 s a) through interference between the graphs of Fig. 10 a and Fig. 10 b. The cross 
section exhibits singularities when the gluon becomes soft or collinear to the initial-state quark/antiquark. 
The soft singularities cancel against those from the virtual photon/gluon contributions in qq annihilation, 
when added along Eq. (9), while remaining collinear singularities have to be absorbed in the PDFs by 
factorization. IR and collinear singularities can be treated by mass regularization. 

For applying the phase space slicing method, the eikonal current has to be modified due to colour 
correlations after the emission of the soft gluon (see Refs. [48-50] and Appendix C for details). Colour 
correlation has to be taken into account also when using the Dipole Subtraction Method; we modified 
the formulae of Ref. [45] accordingly, following the guidelines of Ref. [48]. In Fig. 12 we illustrate the 
comparison between the two methods also for gluon radiation, with good numerical agreement. 
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3.6 q(q) g fusion 

A last class of partonic processes at the considered order is given by (anti-)quark-gluon fusion, 

q(pi) g(p2)^ Q a (ki) Q a *(k 2 ) q (k 3 ), 

q(pi) 9iP2) - Q a (h) Q a *(k 2 ) q(h) . (16) 

This IR finite class contributes at 0(a 2 a) through the interference between the diagrams of Fig. 11 a 
and Fig. lib. Mass singularities arise when the incoming gluon and outgoing (anti-)quark are collinear. 
These collinear divergences are again absorbed into the PDFs. Their extraction has been performed using 
the two methods described above in section 3.2. Explicit expressions for the cross section in the collinear 
region can be found in Appendix C. The actual expression of the subtraction function used in the Dipolc 
Subtraction method is obtained from the formulae in Rcf. [51]. Since those formulae are given there for 
the case of photon-quark splitting we have to consistently redo the color algebra. In Fig. 12 we show the 
agreement between the two methods for the example ug — > u L u L *u. 

In specific cases of SUSY parameters, when kinematically allowed, the internal-state gauginos can be 
on-shell. The poles are regularized by introducing the width of the corresponding gluino, neutralino, or 
chargino. Potential problems related to gauge invariance [52] do not occur here. 



3.7 Factorization of initial-state collinear singularities 

The 0(a 2 a) corrections to partonic cross sections contain universal initial-state collinear singularities that 
can be absorbed into the PDFs choosing a factorization scheme where singularities of relative order 0(a), 
the lowest order PDF fi{x) for parton i (= q, q) is related to the experimentally accessible distribution 
fi(x,fi F ) via [53,54] 



fi(x) = fi(x,n F )\l + 

,1-5. 



aef + a s C F 



-1 + lnS s +\n 2 S s 



In 5 S 



?lh, 



V/4 



+ t— HSs) 

4 7T 



f dz f fx \ I .aef 
+ I -f,{-,P F )^(z) — 



C F a s 



2tt 



2tt 



In 



m?(l - zf 



(17) 



with the factorization scale ii f , C F = $,T F = \ and the electric charge ej. The splitting functions Pa, 
Pig are defined in the usual way, 

P*{z) = ±±^- P lg (z) = 7? + (1 - z) 2 , (18) 



and the functions h and H are given by 
2tt 2 

Ws) = 9 + — + 31nS s - 2 In 2 6 S , H(z) = P u (z) In 



2 1 



+ 2z + 3. 



(19) 



The actual effect of the factorization of the initial collinear singularities is to substitute f%{x) by 
fi(x,fj, F ) m the definition of the quark-antiquark luminosity (2) and thus to obtain a further 0(a 2 a) 
contribution to be added to Eq. (9). This contribution reads: 



E 



1 dT ^hl { r)< 



2ae 



ae 



TO 



dr 



1 < t + ^HSs))do^ 6a6a ,(s) + 



2,0 



2a s C F 



ft j i,i 



q q^Q°-Q 



-1-<5 S 

+ / dz 



+ 



i X 

J Tn 



< K ^ z )-<H{z)\ da 2 L° ( ZS ) + ^K 7 "(z)da 1 L 1 --(zs) 

7T 1 7f w l qq^Q a Q a * x n 1 w qq^Q a Q a ' y 



dr 



dL 



<;.<; 



dr 



(r) 



dL 



<w 



dr 



(r) 



dzPqg(z) 



a s T F 
2tt 



1:i I — f I do 



1,1 



qq^QiQ 



(20) 
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parameter 


SPSla' 


SPS5 


SU1 


SU4 


Af V2 
M 
A 
sign(M) 
tan f3(M z ) 


250 GcV 
70 GeV 
-300 GeV 
+ 
10.37 


300 GcV 
150 GeV 
-1000 GcV 

+ 
5 


350 GcV 
70 GeV 



+ 

10 


160 GcV 
200 GeV 
-400 GeV 
+ 
10 



Table 1: Input parameters in the four benchmark scenarios. 



where xo = (4mj~ a )/s, while K"° !t and k c °"(z) are defined as 

<* = lnS s + In 2 S s + (lnd s + ^ In (^j - 1, (21) 
= PM 

The singularities in k"° u (z) cancel in the sum of the real corrections and of the contribution (20), as 
can be easily checked using the analytic expressions of tier 2 ! 1 - - and tier 2 ! 1 - ~ in the collincar 

J J f qq^Q"Q"j qq^Q<*Q"g 

regions [see Eqs. (39) and (43) of Appendix C]. The remaining singularities of the real corrections are 
exactly cancelled against those in n^ ft and in the virtual corrections. The mass singularities in the last 
line of Eq. (20) are cancelled by those of da 2,1 ~ - and da 2 ' 1 ~ ~ _, as can be inferred from the their 

v ' J qg^Q a Q a *q qg^Q a Q a 'q 1 

analytic expressions in the collinear region [Eq. (45) of Appendix C]. 

For the calculation of hadronic observables we use the MRST2004qed parton distribution func- 
tions [39]. Factorization and renormalization scales are chosen as equal, = fip = itIq a . 

4 Numerical Analysis 

For the numerical evaluation and for illustration of the EW effects, we choose four different benchmark 
scenarios: the point SPSla' suggested by the SPA convention [55], the snowmass point SPS5 [56] char- 
acterized by light stops, and two of the points chosen for detector simulation in the ATLAS "Computing 
System Commissioning" exercise [57]: the point SU1 in the coannihilation region, and the point SU4 
characterized by light SUSY particles. The input parameters Mi/ 2 , M , A , defined at the GUT scale, 
and tan/3 are put together in Table 1. The MSSM input for the actual calculation is obtained with the 
help of the program SPheno [58], together with the program SuSpect [59] as a cross check. The pole 
masses of the squarks of the first generation obtained with the two different codes are shown in Table 2. 
Since the quarks of the first two generations are treated as massless, same-chirality and same-isospin 
squarks are degenerate, therefore we do not show the masses of the squarks belonging to the second 
generation. The difference between the masses provided by the two codes is below 1%. The different 
inputs given by the two codes give rise to a differences in the total cross section of the order of 2 — 3%. 
The standard model parameters are taken from Ref. [60] . 

We introduce the following conventions: 

• We will refer to the sum of (D(a s a), 0(a 2 ) and O{a 2 s o) contributions as "the EW contribution". 

• We will use the quantity 8 to denote the relative EW contribution, defined as S — (Onlo — Clo)/ Clo> 
where O is a generic observable and O nlo is the sum of the LO in Eq. (3) and the EW contribution. 
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SPSla' 


SPS5 


SU1 


SU4 


nl u,R 


548.1 


660.3 


739.7 


412.6 




(545.6) 


(657.4) 


(736.3) 


(411.2) 


mu,L 


565.3 


681.5 


765.6 


420.3 




(562.0) 


(677.5) 


(760.7) 


(418.6) 


m d.R 


547.9 


659.2 


738.0 


413.9 


(545.4) 


(656.9) 


(734.6) 


(412.5) 


m d.L 


570.7 


685.5 


769.6 


427.5 




(567.5) 


(681.8) 


(764.7) 


(425.8) 



Table 2: Pole masses (in GeV) of the squarks of the first generation in the various SUSY scenarios. They 
are obtained using SPheno [58]; those computed with SuSpect [59] are quoted inside the brackets for 
comparison. 

4.1 Different squark species 

Electroweak interactions depend on the hypcrchargc of the squarks, hence the production cross sections 
are flavour and chirality dependent. In this subsection we will study the production of four squark species, 
focusing on the SPSla' point. Since the masses of the light quarks can be neglected, the weak eigenstates 
of the squarks are also the mass eigenstates; thus, in the following, the two squarks of a given flavour are 
distinguished by means of their chiralities, Q a — Q L , Q R . 

Dependence on squark flavour and chirality 

In Tab. 3 we show the integrated hadronic cross section for the diagonal pair production of u L , u R , d L 
and c L . In the case of the production of the squarks of the first generation there is a cancellation beetween 
0{a s a) and 0(a 2 ) contributions. The overall 0(a s a + a 2 ) correction is negative and of the same order 
of magnitude as the 0(a 2 a) one. Since they have the same sign their effect is enhanced. In the case of 
c L production the situation is different: O(a s o), 0(a 2 ), and 0(a 2 a) corrections are positive, 0(a 2 a) 
contribution being the most important ones (see also the discussion below). 

As a general remark, the EW effects are always larger for left-handed squarks. For a given chirality 
and generation, the EW contributions are more important in in the case of up-type squarks. For com- 
parison we also estimate the corresponding NLO QCD corrections using the code PR0SPIN0 [25]; they 
are positive, weakly dependent on the flavour of the produced squarks, and of the order of 47 — 48%. 

Fig. 13 shows the relative EW contribution (right part) in the "cumulative invariant mass distribu- 
tion" er(M inv ), that is the cross section integrated up to the value M inv of the squark-antisquark invariant 
mass. A common feature is that in the low invariant mass region the NLO EW contribution is positive, 
rather steeply decreasing as the invariant mass increases, reaching the plateau at M inv > 2000 GeV which 
corresponds to the total cross section. The left part of Fig. 13 shows the relative size of the individual 
contributions arising from the various channels. The contribution from the gluon fusion channel is always 
positive and dominates at lower values of M inv , wheras the qq annihilation channel part is negative. 
Looking at the relative contributions of the different channels in the high invariant mass region, which 
corresponds to the total cross section, one can understand the origin of the different behaviour of the 
NLO EW corrections in the case of u L u L * and c L c L * production. For up-squark pairs, the 0(a s a) 
and 0{a 2 a) terms are dominated by the the qq annihilation channels, which yield a negative contribu- 
tion; for charm-squark production, however, the O{a s o) [0(a 2 a)\ corrections are dominated by the qj 
fusion [gg fusion] channel and thus positive This shows the key role played by the partonic processes 
Q Q, Q' Q' — ► Q a Q a * , where Q and Q' belong to the same isospin doublet. Indeed, in the case of u L u L * 
production their contribution is negative and the largest out of the qq annihilation channels. In c L c L * 
production they are suppressed by the PDFs of the charm and strange quarks and hence the contributions 
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u R u R * 


u L u L * 


d L d L * 


~L~L* 


0{al) 
O(oe s a) 

0(a 2 ) 
0(a 2 s a) 

8{%) 


(36.83 ±0.03) • 10~ 2 
(-9.00 ±0.01) • 10" 3 
( 2.42 ±0.01) • 10~ 3 
(-3.09 ±0.05) • 10~ 3 
-2.6 


(31.31 ±0.01) • lO- 2 
(-3.54 ±0.01) • 10~ 2 
( 2.39 ±0.01) • 10~ 2 
(-1.05 ±0.01) • 10~ 2 
-7.0 


(25.89 ±0.01) • 10~ 2 
(-3.83 ±0.01) • 10~ 2 
( 3.20 ±0.01) • 10~ 2 
(-7.82 ±0.07) • 10~ 3 
-5.5 


(22.65 ±0.01) • 10~ 2 
( 2.82 ±0.01) • 10~ 3 
( 2.11 ±0.01) • 10~ 3 
( 5.89 ±0.01) • 10~ 3 
4.8 



Table 3: Total cross section for the diagonal pair production of different squark species in the SPSla' 
scenario. Beside the LO contribution, of 0(a 2 ), we show the yields of the different orders contributing 
to the NLO EW corrections. Cross sections are given in pb. 8 is defined according to Sec. 4. 



from the qq annihilation channels are negligible. As a result the overall contribution to total cross section 
is negative at the level of 5% for the left-handed up-squarks, while for the left-handed charm-squarks it 
is of the same order of magnitude but positive. 

The contribution of the 177 channel is independent on the squark chirality, determined only by the 
electric charge of the produced squarks, which makes the gj channel contribution for up-squark pair 
production four times bigger than that for down-squarks. Owing to the mass degeneracy between same- 
chirality and same-isospin squarks the gg fusion channel is independent on the generation of the produced 
squark. 

The invariant mass distribution itself is displayed in Fig. 14 for the various squark species, showing 
also the breakdown into the individual channels. For each squark species, the EW contributions are 
positive in the low invariant mass region and become negative for larger values of M inv , reaching the level 
of 15% for u L squarks. 

Fig. 15 contains the transverse momentum distribution of the squarks. Again, the EW effects are 
more pronounced for left-handed chirality yielding more than 30% negative contributions for large pr- 
As new feature, the LO EW contribution can be positive for low px, especially for the d L case, origi- 
nating from the PDF-enhanced parton process uu — > d L d L * through ^-channel chargino exchange. This 
positive part is practically compensated by the NLO O(a 2 s o) contributions in the qq annihilation channel. 

Dependence on the squark masses 

To study the dependence of the NLO EW contributions on the mass of the squarks, we vary m^jj, setting 
m d r — m u,R and rriu t R = ^,1,(1 + e) with e = 0.03, which is the value at the SPSla' point. The values 
are also taken for the other generations as well as for the sleptons. The other parameters are kept as in 
SPSla'. Each parameter point was checked to satisfy the bounds on SUSY particles from LEP [61,62] 
and Tevatron [63], and the bound on the mass of the light Higgs boson h°, which has been computed 
using FeynHiggs 2.5.1 [64-66]. Moreover, each point fullfills the condition |Ap| < 0.025, where Ap is 
the dominant squark contribution to the electroweak p parameter. 

The relative EW contributions are shown in Fig. 16 for the total cross section, for each of the various 
squark types. The quantity £ displayed in the right panel is the fraction of each the gg fusion and the 
qq annihilation channel in the total cross section, at leading order 0{a 2 s ). The qq channel becomes more 
and more important as mj^ increases. This feature, already pointed out in Ref. [19], is a consequence 
of the i-channel gluino exchange diagrams. The increasing importance of qq annihilation allows a better 
understanding of the particular role of the NLO corrections to the qq channel with increasing squark 
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SPS5 


SU1 


SU4 


0{al) 


(10.62 ±0.01) 


io- 2 


(51.77 ±0.02) 


10~ 3 


(16.14 ±0.01) • 




0(a s a) 


(-1.37 ±0.01) 


• io- 2 


(-7.22 ±0.01) 


• 10" 3 


(-1.45 ±0.01) 


io- 1 


0(a 2 ) 


( 9.11 ±0.01) • 


10~ 3 


( 4.73 ±0.01) • 


10~ 3 


( 10.16 ±0.01) 


io- 2 


0(a 2 s a) 


(-4.83 ±0.03) 


• io- 3 


(-2.75 ±0.02) 


• io- 3 


(-2.61 ±0.01) 


io- 2 


S(%) 


-8.9 




-10.1 




-4.3 





Tabic 4: Same as Tab. 3 but focusing on u L u L * production in different SUSY scenarios. 



masses. Especially for left-handed up- and down-squarks, the NLO EW contributions become more 
important than the LO ones, with effects of more than 20%. In the charm-squark production case qq 
channel is subleading with respect to the gg and gj fusion channels due to the aforementioned suppression 
of charm and strange PDFs. The total sum of the EW contributions is shown in the right panel of Fig. 16. 
For illustration, wc also give an estimate of the formal statistical uncertainty <5 stat = (L (t nlo ) -1 / 2 , 
assuming a luminosity L = 100 fb _1 . 

4.2 Different SUSY scenarios 

Here we discuss the electroweak effects in the different SUSY scenarios mentioned above. As a concrete 
example, we consider the production of u L squarks, with the corresponding masses listed in Table 2. 

In Table 4 we show the total cross section for the aforementioned production process. The LO con- 
tribution and the different orders entering the NLO EW corrections are shown separately. As one can 
see the absolute value of the different contributions decreases as the mass of m,u,L increases, while the 
relative yield of the NLO EW corrections increases with the mass of the produced squarks. In the case 
of the SU1 scenario NLO EW corrections are negative and of the order of 10%. The corresponding NLO 
QCD corrections are estimated using the code PR0SPIN0 [25]; they are of the order of 45 — 50%. 

Fig. 17 contains the cumulative invariant mass distribution, again with the individual and the total 
EW contributions, which show a similar behaviour for all the chosen scenarios. Also the differential 
invariant mass distribution, displayed in Fig. 18, has similar qualitative features in all scenarios. At low 
values, the gluon fusion part dominates and renders the total EW contribution positive. At larger values, 
the contributions from qq annihilation turn the EW contribution to the negative region; thereby the NLO 
part is always of about the same size as the LO part. 

In Fig. 19 we show the transverse momentum distribution in the various cases. Again, their shape 
depends only weakly on the scenario. 

This general situation is only slightly changed when kinematical cuts are imposed, as we find from 
repeating our analysis for an exemplary set of cuts on the transverse momentum and on the rapidity of 
the two squarks, 

p T > 150 GeV, \y\ < 2.5. 

The cut on the rapidity is not effective because the NLO EW contributions to the rapidity distribution 
are very small for \y\ > 2.5. More important is the cut on the transverse momentum. It excludes the 
kinematical region where the largest part of the gluon channel contribution comes from. Moreover, this 
cut suppresses also the contribution of the gj channel and enhances the influence of the qq channel by 
excluding the region with a positive px distribution. As a result, the negative EW contribution to the 
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total cross section is larger than without cuts, as one can see from Fig. 20. 
Dependence on the gluino mass 

Finally we study the dependence of the EW contribution as a function of the mass of the gluino m g , 
with the other parameters kept fixed according to the SPSla' point. Again, the parameter range is in 
accordance with the phenomenological constraints described in the previous subsection 4.1. At LO, the 
gluon fusion channel does not depend on the gluino mass, while the qq annihilation channel contribution 
decreases with increasing m g , as displayed in Fig. 21. In the low m g region the two channel contribute 
equally to the production cross section, while gluon fusion becomes dominant as the mass of the gluino 
increases. The relative EW contributions from the various channels are flat, adding up to a total EW 
contribution from —7 to — 3% for gluino masses between 500 and 2000 GeV. Thereby, in qq annihilation, 
both the tree-level term 0(a s a + a 2 ) and the NLO corrections 0(a 2 a), are practically of the same size. 

5 Conclusions 

We have computed the 0{a 2 a) NLO electroweak contributions to the production of flavour-diagonal 
squark-anti-squark pairs in proton-proton collisions, in combination with the electroweak LO tree-level 
contributions of 0(a s a + a 2 ). 

We have performed an explicit study of the electroweak contributions for each case of the four squark 
species in the first SU(2) doublet, with a numerical analysis for the LHC. The electroweak effects can give 
rise to sizeable modifications in cross sections and distributions, in particular for left-handed squarks. 
Thereby, the NLO terms are significant and have to be considered together with the tree-level contribu- 
tions. They show a strong dependence on the squark masses, increasing their relative influence with the 
mass of the squarks. 

Moreover, we have investigated several SUSY benchmark scenarios and found that the behaviour of 
the electroweak contributions is only weakly dependent on the scenario. Also the gluino-mass dependence 
is weak. In summary, the electroweak contributions in squark-pair production can reach 20-25% in size 
and are thus significant; about half is carried by the NLO contributions. As a final remark we would like 
to mention that the NNLO QCD contributions, of 0(af), can be expected to be of similar size. 
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Appendix 



A Feynman diagrams 

In this Appendix generic diagrams for the various contributions to the different channels are shown. We 
choose the up-squark case as a specific example. In the following we will use the label S° (S) to denote 
all the neutral (charged) Higgs bosons, while V° = 7, Z. 
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Figure 1: Tree-level QCD diagrams for qq — > u a u a * and for 175 — » u a w a 
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Figure 2: Tree-level EW diagrams for — > u a u a * . 
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Figure 3: Lowest-order diagrams for photon-gluon fusion jg — > u a u a * 
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Figure 4: Tree-level diagrams for real photon emission in gg — ► u a u a *-f. 
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Figure 5: One- loop EW diagrams for qq — > ?i a u a *. The diagrams with counter terms can be computed 
according to the Feynman rules in appendix B. The renormalization constants in the counter terms have 
to be evaluated at 0(a). 




Figure 6: One- loop EW diagrams that enter only in the case of the process dd — ► u a u a * . The diagrams 
containing the counter terms can be computed according to the Feynman rules in Appendix B. The 
renormalization constants in the counter terms have to be evaluated at 0(a s ). 



14 




u u a u u a u u a u u a u g u a u Xi u a u Xi u a 




u u u u 



Figure 7: One- loop EW diagrams that enter only in the case of the process uu — ► u a u a * . The diagrams 
in the last row contain the counter terms listed in Appendix B. The renormalization constants in the 
quark-squark-gluino counter term have to be evaluated at 0(a), the other ones at 0(a s ). 
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Figure 8: One loop QCD diagrams for the process qq — ► u a u a * . These diagrams interfere with those of 
Fig. 2 yielding 0(a s a) contributions. The diagrams containing counter terms can be computed according 
to the Feynman rules listed in Appendix B. The renormalization constants appearing in the counter terms 
have to be evaluated at 0(a s ). 



1G 






(b) 



Figure 10: Diagrams for gluon bremsstrahlung from QCD (a) and EW (b) Born diagrams. They con- 
tribute at 0(a 2 s a) through QCD-EW interference 
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Figure 11: QCD (a) and EW (b) Born diagrams for quark gluon fusion channels. Their interference 
contributes at 0{a 2 8 a). 
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B Counter terms and renormalization constants 



Here we list the counter terms for renormalization of vertices and propagators in the one-loop amplitudes 
for squark-pair production. For squarks of the first two generations, we can neglect L-R mixing, and 
weak eigenstates are also mass eigenstates that can be distinguished by their chiralities a = L,R. The 
Feynman rules for the counter terms can be expressed in terms of the field renormalization constants of 
quarks, squarks, gluons, and gluinos, defined from the relation between bare and renormalized fields, 



o b r = g t ( x + \ 6z g) , *r° = *r ( x + \ sz s) > w 

together with the renormalization constants for the strong coupling g s , for the strong Yukawa coupling 
g s , and for the squark masses, which are defined according to 

g^=g T {l + 5Z g ), g b /"=gr(l + 6Z s ), m 2 ^" = m 2 ^ + Sm\ a . (23) 

The actual expressions of the counterterms that are relevant for our squark-pair production processes are 
given below. 

• Vertex counter terms involving gauge bosons: 



JUUUAJUx;'"' = -ig s (5Z Q:a + s -^ + 8Z g )T c (k + k>)^ 



Q" . 
0° 



i9 2 s SZ Q ~J^ + f c ^ A T A )g, v 




-ig s [{^f- + SZ g + 5Z qL ) ltl Lj- + (iffl + SZ g + 5Z qR ) lfl uj + ]T c 




-ie[CY(q)6Z qLl ^- + CX(q)5Z qR ^uj + ] V° = 7, Z 



-ie[CY(Q)5 aL + CX(Q)SaR}8Z Q ~Jk + k')^ V°=j,Z 



k and k' are the momenta of the squark and the antisquark, and they are fixed according to the 
arrow. T c are the color matrices and f ABC the structure constants of the color group. We omit 
the color indices of fermions and sfermions. Moreover, we define 

Cl(q) = e q , C z (q) = ^—(l 3 q -e q s 2 w ), C*(q) = - S -^e q , (24) 



19 



where sw and cw are sine and cosine of the electroweak mixing angle 6w 
Self energy counter terms: 



Q" 

• — x--*- 

Q" 



^ - m l:J sz Q,a - Sm 2 Q - a ] 



iW-m s )5Z s -8mg] 



• Vertex counter terms involving gauginos: 



Q 



( SZg.a + 2SZ d + 5Z g + 8Z QL ) 5 aL 



Q a ( 6Z Q - a + 25 Z g + 6Z S + SZ QR ) 5 aR lo+ } T c 



Q 



U>- 



i%{ 5Z Q - a + 25 Z g + 5Z g + 5Z QR ) 5 aR 
Q ( 5Z Q - a + 25 Z g + 5Z g + 5Z QL ) 5 aL tu+ } T c 



Q 



Q 



0" 



Q" " 



ie[A_(Q) (5Z Q a + 5Z QL )5 aL w_ + A+(Q) (5Z Q a + 5Z QR )5 aR u+] 



ie[A* + {Q) {5Z Q a + 5Z QR )S aR w- + A*_{Q) {5Z Q a + 5Z QL )5 aL 



-i^(SZQ ta + 5Z Q , L )5 aL ^ 



Q' 



Q" 



Q' 



Q° * 



-ie^(SZQ !a + 5Z Q , L )5 aL u + 



The Feynman rules involving Majorana particles follow the prescription of Ref. [67]; in particular 
the fermion flow is fixed according to the arrow depicted in the quark line. As usual, Q' denotes 
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the SU(2) partner of Q. The vertices involving neutralinos contain the quantities 

where Ny is the mixing matrix of the neutralinos. B{Q') can be expressed in terms of the mixing 
matrices U and V of the chargino sector: B(Q') = U* x [B(Q') = V*^\ for up [down] type quarks. 

The renormalization constants of the squark sector are fixed by on-shell conditions (see also Ref. [34,44]), 

{<9S (v 2 ) 1 
' * m i,. = ^{ s a..( ro l,.)}' 
) \p 2 =m 2 - 

k ' / n 2 =m 2 

r Q',a 

where (Q, Q') is either of the two 5J7 (2) doublets (u, d), (c, s), and Eg is the self energy of the squark 

Q a . Due to SU(2) invariance the mass counter term of the left-handed down-type squark is a dependent 
quantity, 

5m %',L = 5m %.L ~ C2 /3 SM w + 4Mw c % s l3 St/3 , (27) 

(where eg = cos#, sg = sin 9 etc. for abbreviation). The counter term 5tp for tan/3 is fixed in the DR 
scheme and can be written in the following way [68,69], 

StfS=U^Kt{Z d AZ z (rn 2 A0 )} , (28) 

where E div denotes the divergent part of the A°Z self energy in dimensional reduction. As pointed out 
in [70], this process-independent condition is also gauge invariant. Furthermore, the W mass counter 
term appears in (28), in the on-shell scheme given by 

8M* v = VKt{Y% v {M* v )}, (29) 

where is the transverse part of the W self energy. 

The field renormalization constants of the quarks are obtained via on-shell conditons as follows [71], 
5Z qa = -<Ke {£,„(m*)} - mpKt (Z qL (p 2 ) + T, qR (p 2 ) + 2T, qS {p 2 )) } (a = L, R) (30) 

with the scalar coefficients in the Lorentz decomposition of the self energy, 

Z q (p 2 ) =jfuJ-T lqL {p 2 ) +jfuj + Y, qR {p 2 ) + m g T, qS {p 2 ). (31) 
Also in the gluino sector we determine the renormalization constants by on-shell conditions, 

5m- g = iwe{m 5 (E Si (m|) + E Sfl (m3) + 2S 5 s(m?))} (32) 
5Z- g = -9le{E^(m?)}-™|lRe{^(E^(p 2 ) + E^(p 2 ) + 2S^(p 2 ))j | 

{Op J 2 =m 2 

K 9 

The renormalization of the strong coupling deserves some particular care. As mentioned in Section 3.3 
the strong coupling g s is renormalized in the MS scheme decoupling the heavy particles (top, gluino and 
squarks) from its runnning. Accordingly, the renormalization constant for g s in (23) is given by [24] 



° Z g = ~t~ 

y 47T 



2 J £-12 \ M 2 j 3 ^2 



(33) 
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where A = 2/e — 7# + ln(47r). The treatment of UV divergences in dimensional regularization violates 
supcrsymmetry at the one-loop level, introducing a mismatch between the strong Yukawa coupling and 
g s . In order to restore supersymmetry in physical amplitudes, cancellation of this extra term is required, 
which at one-loop order can be achieved by modifying the renormalization constant for g s to be different 
from SZ g : 

SZ g = SZ g + ^-. (34) 

For completeness we quote also the field renormalization constant of the gluon in Eq. (22), 

SZ G = 25Z g . (35) 

At 0{a 2 s a) it enters only the one-loop amplitude -Mg^gaQa. > si ncc the gluon only appears in internal 
lines, SZq is cancelled in the sum of self energy and vertex counter terms. 



C Bremsstrahlung integrals 

Here we list the IR and collinearly singular integrals that appear in the phase space integration of the 
bremsstrahlung processes, with either photons or gluons radiated. In the phase space slicing method, 
cuts are imposed: AE = 2S s ^/s on the energy of the emitted photon (gluon), and an angle cut S c on the 
angle between the photon/gluon and the radiating quark via cos 8 > 1 — 5 C . The phase space is thus split 
into a soft and a collinear region that are singular and a complementary non-singular region, which is 
integrated numerically. 

The integration over the soft region can be performed analytically, regularizing the singularities by 
small masses for the photon (gluon) and the light quarks. With the help of explicit formulae for the IR 
integrals [71,72], one obtains the factorized expressions for the cross section given below. In the collinear 
region, the integration can be expressed as a convolution of the lowest-order cross section and a radiator 
function. 



Process gg — > Q a Q a *j 



This process is affected by IR singularities only. Integrated over the soft region, the differential cross 
section reads as follows, 



da 



2,0 



AaAa, = e 2 =(5 F - 5 FF ) da A 



where 



5f = 
S F f = 



In 



/4A£ 2 \ 1 fl-(3 



2 ( S - 2m l,a 
{ 



1 + {4AE 



In 



In 



1-(3J \ A 2 



- Li 2 



2/3 \ _ l ln 2 A + 



1 + 4 \l-f3 



(36) 



A is the infinitesimal mass regularizing the IR divergencies, and (3 = ^Jl — (4m^ J s). 
Process qq — > Q a Q a *j 

The differential cross section integrated over the soft region can be expressed in terms of the 0(a 2 s ) cross 
section for qq — > Q a Q a * and a IR-singular factor, 



da 



2.1 



a 



e 2 q (Si - 5n) + e2 q( s F - Sff) + e q eQ{5 IF - S F i) 
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where e q = | if q = it, c and e q 



\ otherwise. Furthermore, 



(5/ = In 

Sn = In 

(5/ir = In 

- 2 



(4AE 2 



In 



In 



/4A£ 2 



Li 2 1 



In 



7T 

y 

/4A£ 2 

(1 - 0)st 
2m% {m% -t) 

Q,a Q.a 



■In 2 
1 

f 2 



2 ; ' 



In" 



Lio 1 



1-/3 
.1 + /3 

(l + /?)at 
2m 2 - (m 2 - -i) 

Q,a x Q,a 



7T / 

y + 2Li 2 1 



(38) 



si 



5fj can be obtained from #j,f by the substitution t — > u. 

The differential cross section integrated over the collinear region can be written in terms of a convo- 
lution integral, 



da 2 !} - - (s) = — ^ / dz 

qq~*Q a Q a *t^ > 2 7T L 



hi 



( SS ° 

\2ml 



- 1 



P qq {z) + {l-z)\da^ aQa ,{zs), (39) 



with a; = (4m^ q )/s and the quark splitting function P qq {z) from Eq. (18). 



Process qq — > Q a Q a *g 

The singularities affecting this radiative process are Abelian-like, similar to the case of photon radiation, 
and thus can be treated by mass regularization as well. The differential cross section integrated over the 
soft region can also be expressed in terms of 0(a s a) contributions to the cross section for qq — > Q a Q a *, 
but only together with a rearrangement of the color structure. The emission of a gluon as a colored 
particle leads to color correlations in the eikonal current, which can be taken into account following the 
prescription of Ref. [48] , yielding the result 



T 2,1 

qq^Q a Q a 'g 



da 



1,1 

qq^Q a Q" 



5fi — SiFlda '_ 



qq^Q a Q a ' 



+ S n + 6 



da 1 - 

qq->Q a Q a 



(40) 



with Cf — | and N = 3. In order to specify the color-modified "cross sections" da and da, we first 
separate the tree-level amplitudes for qq — ► Q a Q a * into color factors and reduced matrix elements, 
according to the s- and i-channel diagrams in Fig. 1: 



i/O, qcd [ci,C2,C3,C4] 
qq^Q-Q" 



EfrpC rpC yW0,qcd (s) . rjiC rpC lv0,qcd(t) \ 

^ c 2 ci c 3 c 4 y l qq ^QaQa* c 3 Ci c 2 c 4 y qq^Q-.Qa.» ) 



M 



0,cw [Ci,C2,C 3 ,C 4 ] _ 



qq^Q a Q" 



= s, 



0,cw (s) 



+ 4 



0,ow (t) 



(41) 



where T c are the color matrices in the fundamental representation. With this notation we can write for 
the color-rearranged contributions, 



da 



1,1 

qq^QoQ* 



dt 1 



J-1,1 



167TS 2 iV 2 

dt 



2 me 



If A/f 0,qcd [ Cl > C2 > Cl > C3 A A// '™ [ C4 ' C2 ' C4 > C3 ] \ 

X V qq^Q a Q" ) qq->Q a Q" J ' 



167TS 2 TV 2 
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where color summation has to be performed over each pair of equal indices. On top, average over the 
initial helicities is assumed. Owing to the particular color structure, da is different from zero only if q = Q. 

The differential cross section integrated over the collinear region can be written in terms of a convo- 
lution integral similar to Eq. (39), 



d 2,1 (s) — a>s< ~' F 



/ dz 

J Xn 



2.0 



In 



{2m* J 



P qq (z) + (l-z)}da^ aQa ,(zs), (43) 



with da 1 ! 1 ~ ~ instead of da*_ 

qq^Q a Q a ' qq^Q a Q" 



Processes qg -> Q a Q a *q and qg -> Q a Q a *q 

These processes exhibit singularities when the final (anti-)quark is cmmitted off the gluon in the collinear 
region. In that region the differential cross section can be written, in analogy to [54], as follows, 



do 1 ' 1 - - (s) 

qg^Q a Q a *q x 



da. 



2.1 



2ir 

OsTf_ 
2ir 



dz < In 



s(l - z) 2 8 c 



2m^ 



P qg (z) + 2z(l - z) \ da'' 1 ,- (zs) (44) 



qq^QiQ" 



dz < In 



sQ. - z)H c 

2m| 



P qg (z) + 2z(l -z)\da 



1,1 

qq^Q°Q° 



.(*«) 



with the splitting function P qg from Eq. (18). 
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Figure 12: Lowest order partonic cross sections for the process gg — > u L u L *^f (first panel), uu — + u L u L *j 
(second panel), uu — ► u L u L *g (third panel) and ug — ► u L u L *u (fourth panel), computed with the two dif- 
ferent methods. A is defined as A = cr siicing — cr Dipolc . The error bars represent the integration uncertainty. 
The SUSY parameters are those of the SPSla' point [55]. 
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Figure 13: Cumulative invariant mass distribution for different species of squark pairs, denned as the cross 
section integrated up to M inv of the invariant mass of the squark-antisquark pair. The left panels show 
the relative contributions from the various channels, the right ones show the complete EW contribution. 
The SUSY parameter point corresponds to SPSla'. 
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Figure 14: Invariant mass distribution for different species of squark pairs, for the SUSY parameter point 
corresponding to SPSla'. The left panels show the relative contributions from the various channels, the 
right ones show the complete EW contribution. 
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Figure 15: Transverse momentum distribution for different species of squark pairs. Notations and input 
parameters as in 14. 
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Figure 16: Squark-mass dependence of the EW contributions. Total EW contribution (left), individual 
contributions from the various channels (central). The panels in the right column show the relative yield 
of the two channels that contribute at LO. 
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Figure 17: Cumulative invariant mass distribution for PP 
Notations as in Fig. 13. 



"X in different SUSY scenarios. 
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Figure 18: Invariant mass distribution for PP — > u L u L *X in different SUSY scenarios. Notations as in 
Fig. 14. 
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Figure 19: Transverse momentum distribution of the process PP — > u L u L *X in different SUSY scenarios. 
Notations as in Fig. 15. 
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Figure 20: Same as Fig. 17, but with the kinematical cuts defined in section 4.2. 
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Figure 21: Gluino mass dependence of the total (lower left) and of the individual (lower right) EW 
contributions to the total cross section for PP — > u L u L *X . The upper panel shows the relative yield of 
the two channels that contribute at LO. 
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